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We provide a derivation for the particle number densities on phase space for scalar and fermionic 
fields in terms of Wigner functions. Our expressions satisfy the desired properties: for bosons the 
particle number is positive, for fermions it lies in the interval between zero and one, and both are 
consistent with thermal field theory. As applications we consider the Bunch-Davies vacuum and 
fermionic preheating after inflation. 
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I. INTRODUCTION 

The notion of particles is very intuitive, and at the clas¬ 
sical level, in statistical physics, the dynamics is very suc¬ 
cessfully described by the classical Boltzmann equation 
for particle densities in phase space. In quantum physics 
however, the uncertainty principle seems to prohibit the 
use of phase space densities, and they are replaced by 
their closest analogues, the Wigner functions [3, [4]. Yet, 
strictly speaking they can neither be interpreted as par¬ 
ticle numbers nor as probability distributions on phase 
space, since they may aquire negative values. Attempts 
have been made to define particle number in relativis¬ 
tic quantum kinetic theory [3|, but so far there exists no 
result that would be applicable to general situations. 

In spite of those difficulties, the dynamics of quantum 
fields and particle numbers in the presence of temporally 
varying background fields has been extensively studied 
and is well understood 3,00- The particle number 
operator can be calculated by a Bogolyubov transforma¬ 
tion rotating the Fock space to a new basis, which mixes 
positive and negative frequency solutions. 

In the analysis presented in this paper we show that 
the Wigner function, which we take here as an expecta¬ 
tion value with respect to the ground state of the original 
basis, provides the necessary information about the ro¬ 
tated basis to calculate the particle number produced by 
the coupling to time-dependent external fields. 


II. SCALARS 

A. Scalar kinetic equations 

As the first model case, we consider a massive scalar 
field minimally coupled to gravity, such that in a confor¬ 
mal space-time, with the metric of the form gp V = a 2 rg JlV , 
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the Lagrangean is given by 

s/^gC® = ^a 2 if u {d^){d v ^) - ia 4 m|$ 2 , (1) 

where 77^ = diag[l, —1, —1, —1] is the Minkowski (flat) 
metric, and a = a{rf) is the scale factor. For example, in 
inflation a = — 1 /(Hrj) (77 < 0), while in radiation-matter 
era, a = a r ri + a m r] 2 - Here 77 denotes conformal time, a r 
and a m are constants. 

We quantize the theory CD by promoting to an 
operator, 

$0) = ^ = ^7 X e_lk X (yk(??H + ‘P-kWalk) > 

k 

where V denotes the comoving volume. The mode func¬ 
tions obey the Klein-Gordon equation 

(' d 2 a"/a)ipk = 0, (2) 

where ' = d/di 7, = k 2 + a 2 m|(?7) defines the single 
particle (comoving) energy, and we take for the Wron- 
skian 

Throughout this paper we assume that the modes p >k = 
ipk (k = jkj) are homogeneous, which is justified when 
the mass is varying slowly in space, such that we can ig¬ 
nore its gradients. The field <p = a<f> obeys the canonical 
commutation relation, 

[¥>(x, rj), 0„<p(x', 77)] = iS 3 {x - x'), (4) 

which implies [ak, a' k ,] = <5k.k'- 

The fundamental quantity of quantum kinetic theory is 
the two-point Wightman function, which we here write 
for the ground state |0) annihilated by ak, dk|0) = 0. 
With the rescaling suitable for conformal space-times, it 
reads 

iG < {u,v) = a{u)iG < {u,v)a(v) = (0|<p(v)^(zi)10), (5) 

and its Wigner transform is defined as 

iG < (k,x) = j d A re lk ' r iG < {x + r /2,x — r/2), 
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which satisfies the Klein-Gordon equation [3, [§} 

( - ikod v + ^ d 2 - k 2 + mJ(j;)e“^ d ^ dk o^iG< = 0, (6) 

where = a 2 m^ — a"/a. It is then useful to define the 
n-th moments of the Wigner function, 

/n(k, x) = J ^■k'SiG < {k,x). (7) 

Taking the 1st (0th) moment of the imaginary (real) part 
of Eq. © yields 0 

fi - = 0, i/o - h + Wk/o = 0, (8) 

with (Dk = k 2 + to 2 . Eliminating /2 from © yields [§} 


with |a:fe| = y/l + |/3fe| 2 , and where we have introduced 
the phases 


A k = l A kl ex P(^A') > A k = |A_k| exp (i<j> A ) (19) 
ak = \a k \exp(i(/) a ) , f3 k = \Pk\ exp (i<f>p) . (20) 

Eqs. m and (TT51) can be combined to give 


COS(0A + <t>a - M 


(M 2 + |/3 fc | 2 )o k -f4 

2|ttfc| \Pk\ | A k| 


( 21 ) 


while m yields an expression for cos(</>a — 4> a — </>,a). 
Upon squaring Eq. dm and making use of sin 2 ((C) = 
1 — cos 2 (C), we find that 


/o" + 4^k/o + 2(wk) / /o — 0- 
This can be integrated once to give 

rtfZ + lfo"fo-\f? = J, 


(9) 


U 2 -|A k | 2 = f4 2 -|A;j 2 (22) 


is an invariant of the Bogolyubov transformations (1121) . 
(10) Next, we solve (fl8l) for rik = \Pk\ 2 to find 


where the integration constant is obtained by making 
use of fo = |v?k| 2 (cf. Eq. (l27l) below), Eq. © and the 
Wronskian ©• 


n k ± = - 


UkU'k ± y^A^ 


^(^-nj + iAkl 2 ® 2 ) 


2(U 2 - |A k | 2 rc 2 ) 



B. Bogolyubov transformation 


The Hamiltonian density corresponding to the La- 
grangean 0 reads 


H = — ^ |Uk(akQk + Q k a k) + (A k CTkQ-k + h.c.) j 

k 

U k = W k - (a'/a)ip k \ 2 + Wk Wk\ 2 

A k = (^4 — ~+ U; k < / 3 fe' (11) 

Consider now the homogeneous Bogolyubov transfor¬ 
mation 


(2k 

At 


with the norm 


Oik PI 

Pk Oi% 



( 12 ) 


\otk\ 2 — \Pk\ 2 = 1, (13) 

upon which Ak and Uk transform as 

A^ = —2alPk&k + (Ofe) 2 A k + /3 2 Ak (14) 

^k = {\a k \ 2 + \P k \ 2 )fl k -alP* k A k -a k p k At . (15) 

In terms of real and imaginary parts, these equations can 
be recast as 


2|a fe ||/3 fe |U k + |Ak| cos^x + pa-pff) 

— (|ckfe 1 2 +1/3fc | 2 ) | A k | cos(<)>a pet Pp) — 0, (16) 
|Ak| sin(0A'+^a-0/3)- |Ak| sm{(j)K-(j) a -^ =0, (17) 
Hk-(|afc| 2 +l/3fe| 2 )U k 
+2|a fc ||/3 fe ||Ak| cos(0 A -0a-^)=O, (18) 


where x = cos(^>a — <p a — <pp)- Upon extremizing this 
with respect to x , one can show that a maximum is 
formally reached for 2 i 2 lax = Uk/|Ak| 2 , which must be 
greater than one if the Hamiltonian m is to be diago- 
nalizable. Taking account of a; 2 < 1, one finds that the 
maximum for rik± is reached when x 2 = 1, for which 


<WSJ£ - |Aul= + AJJ 2 ± |A k | Ak 1 

” k± - mi - iAkP>-2- (24) 

Since n k - = 0 when | A^. | = |Ak|, the physical branch 
corresponds to rik = rik- ■ Furthermore, when considered 
as a function of |AjJ, 7i k = n k ^ monotonously increases 
as |AjJ decreases, reaching a maximum when |A(J = 0 
(see figure 0 , for which the particle number is 

n k = (Olaj^klO) = - i, (25) 


where o>k = y k 2 + a 2 ?n^. This definition, which cor¬ 
responds to the (constrained) maximum possible particle 
number a detector can observe, we shall use as our defini¬ 
tion for particle number on phase space. Moreover, note 
that, in terms of thus transformed creation and annihi¬ 
lation operators and fik, the Hamiltonian is diagonal 

H = 7777 ^ ^ ^k^k^k A ), [dk; ^k.k' ■ (26) 

k 

such that our definition agrees with the one advocated, 
for example, in Refs. 0 , 1 - 
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FIG. 1: Particle number rik as a function of || for |flk| = 2, 
|Ak| = 1. Provided |A[J < |Ak|, rik maximizes at |A]J = 0. 


This is to be compared with (l2j . which hnds n k ~ 
(— kr/)~ 3 (—< 1). We suspect that tire difference 
is due to the approximate method used in 0- On the 
other hand, when considering the transition from de Sit¬ 
ter inflation to radiation, one finds that the spectrum 
n k ~ (—fcryo) -4 (—kijo) 1 ( rjo denotes conformal time 
at the end of inflation) is produced |l3| . 

As a consistency check, we now apply (1291) to thermal 
equilibrium, where the Wigner function is (c/. Ref. 0 ) 

iG< = 27rsign(fco)6(fc 2 - m l) el 3 ko^ 1 - (32) 

By making use of 0 and (1551) we obtain the standard 
Bose-Einstein distribution, = l/(e^“ k — 1). 

Recently, an interesting particle number definition has 
been proposed in Ref. [151 ], according to which (expand¬ 
ing space-times are not considered): 


C. Particle number in scalar kinetic theory 

It is now a simple matter to calculate the particle num¬ 
ber in terms of Wigner functions. Making use of © 
and 0 we find 

IWfcl 2 = /o, W k ? = \fZ + u>lfo, ( 27 ) 
from which it follows 

= 2 (*2/. + Iyy)-£(£/»). (28) 


(ffk + — |0fc| 2 10/J 2 — fo^fo + w k/o) ■ ( 33 ) 

Note that in adiabatic domain, in which /q' —> 0, 
Eqs. (1551) and (1551) both reduce to rik —> Wk/o — 1/2, such 
that for example in thermal equilibrium of a free scalar 
theory (1551) . both definitions yield the Bose-Einstein dis¬ 
tribution. According to the authors of 0 , the def¬ 
inition tf55D should be applicable to general situations 
(whenever there is a reasonably accurate quasiparticle 
picture of the plasma), and it is obtained as a consistency 
requirement on the energy conservation and quasiparticle 
current relation, respectively, 


We then insert (l28l) into (1551) to get 

This is our main result for scalars, which is positive, sim¬ 
ply because = \/3k\ 2 > 0 (see Eq. (1551) '). Eq. (1551) is of 
course not a unique definition of particle number. Indeed, 
any Bogolyubov transformation 01 corresponds to 
some particle number definition. Our definition (1291) is 
however the special one, in that is correponds to the de¬ 
tector with the best possible resolution, i.e. which mea¬ 
sures the maximum number of particles, as we showed in 
section IIIBI 

We now apply (1551) to the Chernikov-Tagirov [ld| 
(Bunch-Davies [Tit]) vacuum, 

(30) 

which corresponds to the mode functions of a minimally 
coupled massless scalar field in de Sitter inflation, a = 
-l/Hr] (cf. Eq. 0), for which f 0 = (2k)~ 1 [l + l/(kr]) 2 ), 
leading to the particle number 

= 537 = “/ST (31) 


+ \w k \ 2 = wk(^ + ffk), w k |</fc | 2 = \ + ff k - 

(34) 

The consistency is reached when the kinetic and potential 
energies are equal, in which case a generalized quasipar¬ 
ticle energy is given by, = |<(4| 2 /|0 fc | 2 . 

In order to make a nontrivial comparison, consider now 
a pure state of a scalar theory interacting only weakly 
with a classical background field (which can be described 
by a time dependent mass term). The WKB form for the 
mode functions can be recast as 


0 * = ~^= (aoe~ l f V £kW)dri '+ p 0 e l f V 
v2e k 

^ k = -i^( ao e- i r p oe z r ek ^ )dv ')- 


(35) 


le' 

9 —/’fc- 
2 fk 


where e k satisfies, = w k _ (V 2 ) e k/ e k + (3/4)(e"/e k ) 2 . 
In a free theory |a 0 | 2 — \Po\ 2 is conserved, and it is usually 
normalised to one. In an interacting theory however, the 
single particle description breaks down, and consequently 
|Q)o1 2 — \/3o\ 2 is not conserved. For the purpose of this ex¬ 
ample, we assume that the interactions are weak enough, 
such that |ao| 2 — |/?o| 2 is changing sufficiently slow, and 
the subsequent discussion applies. In the adiabatic limit 
£k —> Wk —> constant, the particle number (1551) and (1551) 
of the state (1551) is simply = |/?o| 2 - 
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On the other hand, when applied to the state (l35l) , the 
definition (1331) yields an oscillating particle number even 
in adiabatic regime, 

(fik + |) ~ ^ + (1 + |/3 0 | 2 ) |/?o| 2 sin 2 (2e k 77 - Xa + Xp), 

(38) 

where «o = \ao\e lXa , Po = \Po\e lXf< , which is positive and 
bounded from above by n k < |/3o| 2 = HI]. Hence, 
for the state (13511 our particle number definition (1751) pro¬ 
vides an upper limit for (|33[) . This was to be expected, 
considering that Ecp (1751) was derived in section III Bl by 
an extremization procedure over the Bogolyubov trans¬ 
formations CED- We expect that a similar behaviour per¬ 
tains in other situations. 


III. FERMIONS 

Provided the fields are rescaled as a 3 / 2 ?/) —> if and the 
mass as am —> m, the fermionic Lagrangean reduces to 
the standard Minkowski form, 

ip i<jhp- if(m R + i'y 5 mi)'ip, 

where, for notational simplicity, we omitted the rescaling 
of the fields and absorbed the scale factor in the mass 
term. Note that the complex mass term m = m R (if) + 
imz{rj) may induce CP-violation (c/. Ref. 0 )- 
The fermionic Wigner function, 


where 


foh 

fih 

fih 

fsh 


Tr 

Tr 

Tr 

Tr 


(i p h ) J ^r(-H 0 S<) 

(—/ik • 77 5 P ft ) J ^(—i7°S < ) 

(- ihk-jP h ) J -^-{—i'y°S < ) 
(-7 5Ph ) l 


(41) 


and Ph = (1/2) [1 + hk • 7 y 5 ] denotes the helicity projec¬ 
tor. Eq. (l39ll expresses the conservation of the Noether 
vector current. The traces of the imaginary parts of (1571) 
decouple from (T7U1) at tree level, and hence are of no im¬ 
portance for the analysis presented here. The moments 
f a h can be related to the positive and negative frequency 
mode functions, Uh{k, t) and Vh{k, t) = — iy^ith/k, t))*, 
respectively. They form a basis for the Dirac field, 


^) = ^E e lkx (^akh + v h b ] _^ h ), 

kh 



® 6m 


where 6i is the helicity two-eigenspinor, h£h = h£h- The 
Dirac equation then decomposes into 

id 0 L h - h\k\L h = m R R h + imiR h 
id 0 R h + h\k\R h = m R L h - imiL h . (42) 


iS < (k, x) = — J d 4 re lk ' r (0\if(x — r/2)if(x + r/2)|0) 

satisfies the corresponding Dirac equation which, in the 
Wigner representation, reads 

(ft + ^7°^t - (m R + i'y 5 m I )e~^ tdk «^iS < = 0, (37) 

where (i'y°S < )^ = i"f°S < is hermitean. The helicity op¬ 
erator in the Weyl representation h = k- 7 ° 77 5 commutes 
with the Dirac operator in (EZD, such that we can make 
the helicity block-diagonal ansatz for the Wigner func¬ 
tion (c/. Ref. 0 ) 


Note that these equations incorporate CP-violation and 
thus generalize the analysis of Refs. S 0. Now, 
from 611) one can derive 439D-6QD by multiplying with 
Lh and Rh and employing 

foh = \Lh\~ + \Rh\ 2 , f?,h = \Rh \ 2 ~ \Lh \ 2 
fih = -m{L h R* h ), f 2h = 2%{L* h R h ). (43) 

The Hamiltonian density reads 

k h k/i) “1“ (-A-kft.fr— kft&kftll.C.)j* 

kft 

where 


iS K = ^2 iS h > ^*7o S£ = i(l + hk ■ it) <8 p a g a h, 

h=± 

(38) 

where k = k/|k| and cr a , p a (a = 0,1, 2, 3) are the Pauli 
matrices. Taking the traces of {1, — /i 7 * 7 5 , — ih'f, — 7 5 } 
times the real part of (l37l) , and integrating over k 0 , yields 
the kinetic equations for the 0th momenta of g a h , 

foh = 0 (39) 

fih + 2h\k\f 2 h - 2 nnf 3h = 0 
fih - 2h\k\fi h + 2m R f 3h = 0 
hh - 2m R f 2 h + 2m!fi h = 0 , (40) 


n kh = hk(\L h \ 2 -\R h \ 2 )+mL* h R h +m*L h R* h 
A kh = 2 kL h R h - hm*L 2 h + hmR 2 h , (44) 

with {a k /,, aj^,} = hh,h'&k,k') {^k/M = 6). 7 ' 5'k.k' ■ 

We now use the Bogolyubov transformation 


( a kh \ = f a k h Pkh \ ( akh \ 

U-kJ V-^^JWtkJ 


to diagonalize the Hamiltonian, where a k h and /3 k h are 


1 / akh 

2 \ p kh 


Pkh 

a k h 




\a k h\ 2 + \Pkh\ 2 


1, (45) 
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leading to the particle number density on phase space, 

(46) 


|2 1 h 


n kh - \p kh \ - 2 2wk , 


where now to k = y/k 2 + \m\ 2 . 

To construct the initial mode functions in the adiabatic 
domain, rj —► — oo, we use the positive frequency solution 
and its charge conjugate, 




a o L h + Po L h \ 
a 0 R+ + PoR h ) ’ 


l«ol“ + I/?o I 2 — i- 


From the Dirac equation under adiabatic conditions it 
follows 


Lt 


K 


/u>k + hk 

T ■ m f 

V 2o; k ’ 

h ~ 1 \m\ V 


2 ^'k 


•\/2wk(wk + hk) 


L-h — * 


\/ 2 w k (uk - hk) 


These mode functions correspond to an initial particle 
number = |/? 0 1 2 ■ We now make use of (l43l) to express 
flk/i in terms of the Wigner functions, 

^k h = ~{hkf 3h + m R fi h + m I f 2 h)i (47) 
which implies our main result for fermions, 

n-kh = 7 T— (hkf 3h + rn R fi h + mif 2 h) + 7 : • (48) 

ZCt/k ^ 


Note that in the limit m —> 0, this expression reduces 
to the phase space density of axial particles. Moreover, 
0 < n k h = \Pkh\ 2 = 1 ^ \oikh\ 2 < 1 (see Eqs. (I45H46D ). 

As an application of E q. (14811 we consider particle pro¬ 
duction at preheating [9j, |l7[, in which the fermionic mass 
is generated by an oscillating inflaton condensate. As¬ 
suming that the inflaton oscillates as a cosine function 
results in a fermion production shown in figure [2] Ob¬ 
serve that, even for a relatively small imaginary (pseu¬ 
doscalar) mass term, particle production of the opposite 
helicity states is completely different, implying a nonper- 
turbative enhancement of a CP-violating particle density, 
?ik+ — «k-, which may be of relevance for baryogenesis. 

When applied to thermal equilibrium, where (c/. 
Ref. [3) 

iS < = -Qi +m R - i^m I )5(k 2 - \m\ 2 ) , (49) 

we find 

foh = 1, 

fih = {2m R / w k ) [{exp(/3cc k ) + 1} _1 - 1/2], 
fih = (2m//a) k )[{exp(/3a;k) + 1} _1 - 1/2], 
f 3h = (2hk / Uk)[{exp(pujk) + 1} _1 - 1/2], 

such that Eq. (l48l) yields the Fermi-Dirac distribution, 

nkh = l/(e /3wk + 1). (50) 



FIG. 2: The number of produced fermions as a function of 
time with helicity h = + (solid) and h = — (dotted), mass 
m/uJi = 10 + 15cos(2t) — isin(2r), |k| = uji, t = unt, where 
uii denotes the frequency of the inflaton oscillations. 

IV. MULTIFLAVOUR CASE 


We now generalize the definition of particle number 
in terms of two-point functions to the case of several 
species, mixing through a mass matrix. While in the 
single flavour case always an equal number of particles 
and antiparticles is produced, we will here encounter the 
creation of a charge asymmetry when the mass matrix 
is nonsymmetric. Because of this charge violation, the 
orthogonality of particle modes with respect to antipar¬ 
ticle modes is not preserved under time evolution, and it 
is thus impossible to expand the field operators in terms 
of an orthogonal basis. 

Hence, the use of the basis-independent two-point 
functions is advantageous. We can either calculate the 
time evolution of the system in terms of these quanti¬ 
ties or measure them, since they correspond to physical 
charge- and current densities. When finally the mass 
matrix is diagonal and only adiabatically slowly evolv¬ 
ing, there exists a well-defined basis, in terms of which 
the Hamiltonian is diagonal. We use this basis to de¬ 
fine the particle number operators and construct their 
expectation values out of the two-point functions. 


A. Fermions 


Since Dirac spinors naturally include particle and an¬ 
tiparticle modes, we first discuss here the fermionic case. 
We decompose the mass matrix M into a hermitean and 
an antihermitean part, 

M h = \(M + Aft), M a = i(M - Aft), (51) 

Z ll 

such that the Dirac equation reads 

[i$ - Mh - h b M A ] tj = 0 . 


(52) 
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One can then attempt to proceed as in the single flavour wher a, b are flavour indices. These obey the equation of 
case and to construct the field operators as motion 


,,-ik-: 


H x ) = —y— tr is( k > t Kj( k ) + V h ij ( k . t) b i i (- k ) ( 53 ) 


k hj 


hj 


^ X ) = E V hji M 


khj 


with the mode function 


Uhii = ' Lhi3 


L 

Rh ij 




(54) 


(# + °dt - ( M h + i^M A )e-^ d ^iS< = 0. (60) 

As described for the single flavour case in section Hill this 
can be simplified and yields 

foh + i [Mr, flh] + i [Ma, f2h] = 0 
fih + 2h\k\f 2H + * [M h , f oh ] - {M a , f 3h } = 0 
hh ~ 2h\k\fih + {M H , f3h} + i[M A , foh] = 0 

hh — {M H , f2h} + {M A , fi h } = 0.(61) 


and its charge conjugate 

v hij = -n 2 (u hij y = cu hij c~ l = ( 

(55) 

This procedure however fails, when M is not a sym¬ 
metric matrix, which can easily be seen by plugging Uhij 
into Eq. (1521) 

{id 0 - /i|k|} luLhij = M^Rhij + iM^Rhij 

{ido + ft|k|} laRhij = MgL hlj -iM£L h ij (56) 
and Vhij, respectively, 

{id 0 — /i|k|} t a L hi j = M”*R hlj +iM A *R hlj 
{id o + h\k\} t a R h ij = M%*L hlj - iM$*L hlj ,(57) 

where summation over the repeated index l is implied. 

Obviously, when M is not symmetric, Eqs. AMD 
and (1571) are inconsistent. In particular, for nonsymmet- 
ric M, the orthogonality condition 

Ul il Vsij= o, (58) 

is not preserved at all times, and hence, the expansion 
of the field operators (1551) is not suitable. This complica¬ 
tion can however lead to the generation of a net charge 
stored in the produced particles, because the operation 
of charge conjugation becomes time dependent, an effect 
which may be of relevance for baryogenesis [18l |. 

The construction of an appropriate Bogolyubov trans¬ 
formation for the case of a symmetric mass matrix is 
discussed in Ref. fl9j . In comparison with the single 
flavour case this procedure is fairly complicated. For the 
general case, we therefore refrain from a computation of 
a Bogolyubov transformation and the time evolution of 
Heisenberg creation and annihilation operators. 

It is more convenient to calculate the time evolution 
of the initial state in terms of two point functions. We 
straightforwardly generalize the formalism for the single¬ 
flavour Wigner funtions to the multiflavour case by defin¬ 
ing 

iSfj{k,x) = - J d 4 re lk ' r ( / ipj(x - r/2) r tp i {x + r/2)), (59) 


As already noted in [l8|, we can infer from these equa¬ 
tions as a necessary condition for the nonconservation of 
the charge density foh, that M must not be symmetric, 
in accordance with our discussion above. 

Now assume, that after some time evolution, M has 
become symmetric and slowly varying. Then, it is pos¬ 
sible to expand the field operators as in Eq. (1531) and to 
define the expectation values of the number of particles 


Khi = ( a L( k Hi(k)) 

(62) 

and antiparticles 


Khi = (&L( k )M k )). 

(63) 


Moreover, we choose this basis such that the Hamilton 
operator is diagonal and reads 

H = ^Y. ( ft |k| L\L h + L{[M H + %M A )R h 

k hij 

- h\k\R[R h + R{ [M h - iM A ] L h ) 

/ ij 

x ( a L( k Hj( k ) “ M k )klj( k ))- (64) 

We can now also express the functions /A employing 
this basis. Explicitly, they read 

foi(x, k)=-/ d 4 r e lk ' r (’ip h j(x - r/2)'y°ip hi (x + r/2)) 
=(Lfl{+RfR{) x (<4(k)a w (k) +& w ,(k)&L(k)) 
flh(x, k)=-/ d 4 r e lk r (y hj (x - r/2)ip hi (x + r/2)) 
=-25ft(Ll'i?f*) x (a{ v {k)a hl {k) + b hl }k)b{ l {k)), 
f 2 h( x >*■)=-J d 4 re lk " r {^ h j{x-r/2){-i'y 5 )'ip hi {x + r/2)) 
=23 (if i?f ) x (4,,(k)a w (k) + b hl }k)bl(k)) , 

flh( x ’ k )=-J d4r e lk r {i>hj (X - r/2) r y°'y 5 ip h i{x + r/2)) 

=(ifLf -R%*R%') x (al,(k)a hl (k) + b hl ,(k)b{ l (k)) . 
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By comparing with the expression (1641) . we obtain 


(H) = -y ^2 h\k\f 3hii + Mgf lh>ii + M^f 2hii . (65) 

khi 

We define u>i( k) = (k 2 + {Mu] 2 ) 1 / 2 , and since we as¬ 
sumed diagonality of the Hamiltonian, this has to equal 


( H ) = ^X)^( k )(4i( k Ki( k )- b hi( k ) 6 L) 

k hi 

= ( n ihi+ n khz -!). ( 66 ) 

k hi 

while the charge is 

fohu = (aii(k)a hi (k) + 6 h i(k)6j li > = n+ hi - n^ hi + 1. 

(67) 

We thus find the following generalization of (1151) 


+ _ h\k\f 3 hii + M^flhii + M^f 2 hii 

n * hi ~ &*(kj 

- _ h\k\f 3 hii+MVfihii + M^f 2 hu 

Hkhi ~ 2^(k) 


+ ~fohii (68) 
- 2 /o/iii + l)(69) 


can in general be satisfied only if M 2 is real or, more 
precisely, real symmetric. 

Let us therefore assume again, that we are in a final 
state with diagonal and only nonadiabatically varying M. 
We define 


a(k) 

b( k) 


—p [o 1 (k) + m 2 (k)] and 
v 2 

-j= [ ftl ( k ) -*a 2 (k)] . 


Then, we find the charge operator to be 

Q i (k) = (oJ(k)a i (k)-6j(k)6 i (k)) 


(73) 


(74) 


and the Hamiltonian 

H = ^^fli(k) (ot(k)o*(k) +bl(k)bi(k) + l) , (75) 


where f2»(k) = |^'(k) - (a'/a)<pi( k)| 2 + w, ; (k) 2 |(/?j(k)| 2 , 
w 2 (k) = k 2 + a 2 M%. 

We define the multiflavour Wightman function as 

= (¥>]-(w)<Pi(w)> (76) 


which is of course the anticipated result, since the number 
of particles is just the half of the total particle number 
(particles plus antiparticles) plus half of the total charge 
(particles minus antiparticles). 


B. Scalars 


Consider now a complex scalar field <h,; describing N 
flavours, which we expand into its hermitean and anti- 
hermitean parts as follows, 

= +*$?), (70) 

such that the multiflavour field operator is 


4, = ^ = -L 


Vi 

a 


Ve- ik ' : 

aV ^ 

k 


X (vtj ( k , v)Oj ( k ) + ( k , v)a 2 j ( k ) 

+ <filj (-k, v)a]H - k ) + (-k, r/)a 2f (-k)^ , (71) 


where the rescaled fields obey the generalized Klein- 
Gordon equation 

d 2 + k 2 + M 2 - — \ = 0. (72) 

a J il 

Note that this is independent of whether a = 1 or a = 2, 
which is just as in the fermionic case, where the functions 
U and V both satisfy the Dirac equation. The individual 
components and <f> 2 are imposed to be hermitean. 
Therefore, YljiPiji k ) + Viji— k)) has to be real, which 


and adapt straightforwardly the definition of the mo¬ 
menta from the single flavour case. These then satisfy 
the system of equations 

^/o-/2 + ^{M 2 ,/o}+(k 2 -^/o=0 

tf-l[M 2 ,fo] =0 

fi~\ [M 2 , fi] — ^ {(k 2 + M 2 — a"/a)', / 0 } = 0.(77) 

We find Qi( k) = fm(x,k) + 1, which is also in accor¬ 
dance with the [/(l)-Noether charge. Together with the 
identities m , this leads us to 


n 


+ 

k i 


n 


k i 


f» 

, , f I ■> O ii 

- 

4 UJi 

f" 

, , f I JQii 
UiJOWT— - 

4 UJi 


1 d 

2 uji di] 

1 d 

2 u>i dij 


^ — foii'j + yflii (78) 

(^/Oii) - 1(79) 


where njG is the number of particles, n^ i the number of 
antiparticles, and the same simple interpretion as in the 
fermionic case applies. 


V. DISCUSSION 

We have derived general expressions for the particle 
number densities on phase space for single scalars (]29l) 
and fermions (H51) in terms of the appropriate Wigner 
functions. We have then generalized our analysis to 
the case of mixing scalars ([78] [79]) and fermions m- 
E2). All of these expressions are positive, and more¬ 
over, the number of fermions is bounded from above 




by unity , as required by the Pauli principle. In or¬ 
der to incorporate the effect of self-energy into (12!)1) 
and (l48l) . one needs to include this correction into the 
dispersion relation, u = u?(k,x) —> u> + £#(k, x), where 
S ff (k,x) = f[dk 0 /(27r)](l/2)[Z r (k,x)+E a (k,x)], and £ r 
and £“ denote the retarded and advanced self-energies, 
respectively [20(. When the single particle picture breaks 
down it is not clear whether a sensible definition of 
particle number can be constructed. Our analysis can 
be quite straightforwardly extended to include (time- 
varying) gauge fields by coupling them canonically to 
scalars and fermions. 

The kinetic theory definition of the particle number is 
of course and by construction identical with the definition 
in terms of Bogolyubov transformations. The number of 
individual particles is the total energy of the system di¬ 
vided by the energy of an individual particle. Taking the 
point of view of kinetic theory proves advantageous when 
considering the multiflavour case or statistical systems, 
such as the thermal equilibrium. 

While the fermionic particle number definition (H51) is 
generally applicable, the scalar one (El fails however 
when u;j( = k 2 + a 1 2 3 4 5 6 7 to| < 0, which can happen at phase 
transitions. Then Ok < |Ak| in (fill) , and the Bogolyubov 
transformation m does not have a solution. Neverthe¬ 
less, even in this case, the energy density on phase space 
flk in Eq. (1251) is well defined, and should be considered 
as a fundamental quantity of kinetic theory. Another 
important quantity is A£ = (k, — k|if|0), the transition 
amplitude for particle pair creation with the momenta 
{k, —k}; and likewise Ak is the transition amplitude for 
pair annihilation. The appropriate description in this 
case is in terms of squeezed states. For an account of the 
inverted harmonic oscillator in terms of squeezed states 
see eg. Ref. p2j . 

Our definition of particle number can be used for 


studies of quantum-to-classical transition, decoherence 
and entropy calculations of eg. cosmological perturba¬ 
tions [22, [23, [24|]. Moreover, when suitably normalized, 
the particle density rik can be used to define a density 
matrix on phase space, £>k = rik/ ]Ck' n k'- 

In the derivation of our results, we considered pure 
quantum states, yet showed explicitly their applicabil¬ 
ity to thermal states. More generally, our definitions 
are valid if one requires the density matrix g to satisfy 
(akflk)e = ( a k a k)e = 0- These relations hold e.g. for 
eigenstates of the particle number operator = a^a k, 
and, as pointed out in Ref. [2lj . for random phase states, 
a special case of which is the canonical ensemble. States 
of this kind can be treated as a linear superposition of the 
particle number eigenstates which we considered above. 

Finally, we note that after the first version of this ar¬ 
ticle appeared, an out-of-equilibrium investigation of the 
dynamics of chiral fermions coupled to scalars was stud¬ 
ied in Ref. [25[. In order to show that at late times the 
system thermalizes to the Fermi-Dirac equilibrium, the 
authors used a particle number definition, which can be 
in our notation written as 

Ak = ^ ”kh , n kh = i (1 + hf 3h ) . (80) 

Z h=± 

This definition corresponds to the massless fermion limit, 
m —> 0, of our definition (H51) . 


Acknowledgements 

We would like to acknowledge discussions with Jurgen 
Berges, Dietrich Bodeker, Kirnmo Kainulainen and Julien 
Serreau. 


[1] E. P. Wigner, “Quantum Corrections For Thermody¬ 
namic Equilibrium,” Phys. Rev. 40 (1932) 749. 

[2] M. Hillery, R. F. O’Connell, M. O. Scully and 
E. P. Wigner, “Distribution Functions In Physics: Fun¬ 
damentals,” Phys. Rept. 106 (1984) 121. 

[3] Y. Kluger, E. Mottola and J. M. Eisenberg, “The quan¬ 
tum Vlasov equation and its Markov limit,” Phys. Rev. 
D 58 (1998) 125015 hep-ph/9803372 . 

[4] L. Parker, “Quantized Fields And Particle Creation In 
Expanding Universes. 1,” Phys. Rev. 183 (1969) 1057. 

[5] S. G. Mamaev, V. M. Mostepanenko and V. M. Frolov, 
“Production Of Fermion Pairs By A Nonstationary Grav¬ 
itational Field,” Sov. J. Nucl. Phys. 23 (1976) 592. 

[6] D. J. Chung, E. W. Kolb, A. Riotto and I. I. Tkachev, 
“Probing Planckian physics: Resonant production of 
particles during inflation and features in the primor¬ 
dial power spectrum,” Phys. Rev. D 62 (2000) 043508 
hep-ph/9910437 . 

[7] P. f. Zhuang and U. W. Heinz, “Equal-time hierarchies 
for quantum transport theory,” Phys. Rev. D 57 (1998) 


6525 hep-ph/9801221 . 

[8] D. Bodeker, K. Kainulainen and T. Prokopec (unpub¬ 
lished) . 

[9] L. Kofman, A. D. Linde and A. A. Starobinsky, “Towards 
the theory of reheating after inflation,” Phys. Rev. D 56 
(1997) 3258 hep-ph/9704452;. 

[10] N. A. Chernikov and E. A. Tagirov, “Quantum Theory Of 
Scalar Fields In De Sitter Space-Time,” Annales Poincare 
Phys. Theor. A 9 (1968) 109. 

[11] T. S. Bunch and P. C. Davies, “Quantum Field Theory 
In De Sitter Space: Renormalization By Point Splitting,” 
Proc. Roy. Soc. Lond. A 360 (1978) 117. 

[12] M. Mijic, “Particle production and classical condensates 
in de Sitter space,” Phys. Rev. D 57 (1998) 2138 [gr- 
qc/9801 094]; 

[13] L. H. Ford, “Gravitational Particle Creation And Infla¬ 
tion,” Phys. Rev. D 35 (1987) 2955. 

[14] M. Le Bellac, “Thermal Field Theory,” Cambridge Uni¬ 
versity Press (1996). 

[15] G. Aarts and J. Berges, “Nonequilibrium time evolution 



9 


of the spectral function in quantum field theory,” Phys. 
Rev. D 64 (2001) 105010 hep-ph/0103049 . 

[16] K. Kainulainen, T. Prokopec, M. G. Schmidt and S. We- 
instock, “Semiclassical force for electroweak baryogen- 
esis: Three-dimensional derivation,” Phys. Rev. D 66 
(2002) 043502 hep-ph/0202177 ; 

K. Kainulainen, T. Prokopec, M. G. Schmidt and S. We- 
instock, “First principle derivation of semiclassical force 
for electroweak baryogenesis,” JHEP 0106 (2001) 031 
hep-ph/0105295 . 

[17] M. Peloso and L. Sorbo, “Preheating of massive fermions 
after inflation: Analytical results,” JHEP 0005 (2000) 
016 hep-ph/0003045 . 

[18] B. Garbrecht, T. Prokopec and M. G. Schmidt, “Coher¬ 
ent baryogenesis,” Phys. Rev. Lett. 92 (2004) 061303 
arXiv:hep-ph/0304088 . 

[19] H. P. Nilles, M. Peloso and L. Sorbo, “Coupled fields 
in external background with application to nonther- 
mal production of gravitinos,” JHEP 0104 (2001) 004 
arXiv:hep-th/0103202 . 

[20] T. Prokopec, M. G. Schmidt and S. Weinstock, “Trans¬ 
port equations for chiral fermions to order h-bar and elec¬ 
troweak baryogenesis,” arXiv:hep-ph/0312110 

[21] H. E. Kandrup, “Particle Number And Random Phases,” 
Phys. Rev. D 38 (1988) 1773. 


[22] A. Albrecht, P. Ferreira, M. Joyce and T. Prokopec, “In¬ 
flation and squeezed quantum states,” Phys. Rev. D 50 
(1994) 4807 astro-ph/9303001 . 

[23] D. Polarski and A. A. Starobinsky, “Semiclassicality 
and decoherence of cosmological perturbations,” Class. 
Quant. Grav. 13 (1996) 377 gr-qc/9504030 . 

[24] R. H. Brandenberger, V. Mukhanov and T. Prokopec, 
“Entropy of a classical stochastic field and cosmolog¬ 
ical perturbations,” Phys. Rev. Lett. 69 (1992) 3606 
astro-ph/9206005 . 

[25] J. Berges, S. Borsanyi and J. Serreau, “Thermalization 
of fermionic quantum fields,” Nucl. Phys. B 660 (2003) 
51 arXiv:hep-ph/0212404 ; 

[26] Note that the definition of the quasiparticle energy cj k = 
|<)4|/|0fc|, oscillates even in the adiabatic limit, with 
the minimum and maximum values given by cuk.min = 
ek(|«o| - |/3o|)/(|ao| + |/3o|) and aj k ,m a x = e k /W,min, re¬ 
spectively, such that uju ^ £k in general. This indicates 
that imposing instantaneous equality of the potential and 
kinetic energies may not be appropriate in general situa¬ 
tions. When particle number is understood as an average 
over the characteristic oscillation period however, impos¬ 
ing equality of the potential and kinetic energy may lead 
to a reasonable definition for the particle number density. 




